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Abstract 

Quantum control could be implemented by varying the system Hamiltonian. 
According to adiabatic theorem, a slowly changing Hamiltonian can 
approximately keep the system at the ground state during the evolution if the 
initial state is a ground state. In this paper we consider this process as an 
interpolation between the initial and final Hamiltonians. We use the mean value 
of a single operator to measure the distance between the final state and the ideal 
ground state. This measure resembles the excitation energy or excess work 
performed in thermodynamics, which can be taken as the error of adiabatic 
approximation. We prove that under certain conditions, this error can be 
estimated for an arbitrarily given interpolating function. This error estimation 
could be used as guideline to induce adiabatic evolution. According to our 
calculation, the adiabatic approximation error is not linearly proportional to the 
average speed of the variation of the system Hamiltonian and the inverse of the 
energy gaps in many cases. In particular, we apply this analysis to an example in 
which the applicability of the adiabatic theorem is questionable. 

Keywords: Quantum control; Adiabatic theorem; Interpolation of Hamiltonian 


1 Introduction 

Adiabatic process is aimed at stabilizing a parameter-varying quantum system at 
its eigenstate. This process has many applications in the engineering of quantum 
systems [1, 2, 3, 4, 5], and in particular plays the fundamental role in adiabatic 
quantum computation (AQC) [6, 7, 8]. The adiabatic theorem [9, 10] states that 
a system will undergo adiabatic evolution given that the system parameter varies 
slowly. 

Quantifying the applicability of adiabatic approximations is an interesting topic 
of current research efforts. On the one hand, this kind of research has been spurred 
by so-called shortcuts to adiabaticity [11], and on the other hand recent insights 
from thermodynamics haven put adiabatic processes back into focus [12, 13]. In 
particular, the validity of the adiabatic theorem has been under intensive studies 
both theoretically and experimentally since it was proposed, and much of these 
efforts were devoted to the rigorous description of the sufficient quantitative con¬ 
ditions of adiabatic theorem, and the estimation of the error accumulated over 
a long time [10, 14, 15, 16]. Once the exact knowledge on the adiabatic pro¬ 
cess is available, it is straightforward to apply the results to the optimal design 
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of adiabatic control on specific systems [17, 18]. The most interesting progress is 
that the validity of the adiabatic theorem itself has been challenged in the recent 
decade [19, 20, 21, 22, 23, 24, 25, 26, 27], both by strict analysis and counter¬ 
examples. According to these findings, the errors induced by the adiabatic approxi¬ 
mation could accumulate over time despite certain quantitative condition is satisfied 
[19, 20, 21, 24, 25], e.g., when there exists an additional perturbation or driving that 
is resonant with the system. Particularly as indicated in [24], it is not new that reso¬ 
nant driving can cause population transfer between eigenstates. Also, a proof can be 
found in [25] stating that only a resonant perturbation whose amplitude gradually 
decays to zero can result in a violation of a well-known sufficient condition. 

In this paper we consider the following process: the process starts at t = 0. The 
system Hamiltonian at t = 0 is 77i, and the system Hamiltonian at t = T is 
H 2 = Hi + XAH,X > 0. A is a dimensionless quantity. AH is a fixed operator and 
so the direction of the variation is fixed. We assume Hi, H 2 , and AH are bounded 
operators throughout this paper. T is the evolution time. The transition of the 
system from Hi to H 2 can be described using an interpolating function f{t) so that 

H{t) = Hi + f{t){H2 - Hi) = Hi + Xf{t)AH, (1) 

with /(O) = 0 and /(T) = 1. We work under the condition that a valid perturbative 
analysis of the system evolution is available. This often means A should be smaller 
than a threshold value. It is worth mentioning that the classical adiabatic theorem 
was proved also using a perturbative analysis, which cannot be applied directly to a 
large variation of Hamiltonian. Therefore, our analysis in this paper is not concerned 
with the adiabatic evolution for a large variation of Hamiltonian. However, our 
analysis provides a rigorous estimation of the error accumulated during this small- 
variation evolution for an arbitrarily given interpolation. 

Our work is different from the previous works in two ways. First, instead of study¬ 
ing the evolution of the eigenstates and their corresponding probability amplitudes, 
the mean value of a Hermitian operator is defined as a measure of the error. For 
example, in the context of adiabatic quantum computation where one wants to pre¬ 
pare the ground state of a target Hamiltonian H 2 > 0 whose ground-state energy is 
0, e = {H 2 )pt serves as a good measure of the distance between the real-time state 
Pt and the ground state. This measure resembles the excitation energy or excess 
work performed during the process, as studied in thermodynamics [12]. In this pa¬ 
per we only consider the error accumulated over the entire process, which means 
we are only interested in e = {H 2 )pt- The second difference is that the error, or the 
excitation energy or excess work performed during the process, can be estimated 
with a sufficient precision for arbitrarily given interpolating functions. As a result, 
the parameters which are related to the suppression of the error can be easily identi¬ 
fied. For example, we have e = 0( as A 0 in the case of linear interpolation. 
Here X 2 is the energy gap between the ground and first-excited states of the ini¬ 
tial Hamiltonian. However for the interpolation in the counterexample [19, 25], the 
scaling of e is not so simple. 

This paper is organized as follows. In Section 2, we introduce the model of this 
paper. In Section 3, we give the estimation of the error for linear interpolation. In 
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section 4, we present the general algorithm to estimate the error for an arbitrarily 
given interpolating function. We discuss three examples in Section 5. Conclusion is 
given in section 6. 

2 Definitions and Preliminaries 

The system is defined on an iV-dimensional Hilbert space. We set Dirac constant 
h = 1. II ' II denotes the matrix norm. Two real functions fi{x) and f 2 {x) can be 
denoted as fi{x) = 0 {f 2 {x)), a: ^ oo, if and only if there exists a positive real 
number M and a real number xq such that |/i(a;)| < M|/ 2 (a;)|, x > Xq, where | • | 
denotes the absolute value. 

Let {uji : i = 1,2, ...iV} be the monotonically increasing sequence of eigenvalues 
of Hi, so that Wi > ujj when i > j, and {|*)} be the corresponding eigenstates. 
We denote the energy gap between the ith eigenstate and the ground state as 
Xi = LUi — uJi- Similarly, we define the increasing sequence of eigenvalues of H 2 , 
{Wj : i = 1,2, ...TV} and {Aj}, {|T )} correspondingly. 

For convenience, we also introduce two offset Hamiltonians, Hi and H 2 . The 
Hamiltonian Hi is defined as Hi — Hi — wi, i.e., by offsetting the Hamiltonian of 
the system at f = 0 by a constant operator uji so that Hi > 0. By Hi > 0 we mean 
Hi is positive semidefinite and its the smallest eigenvalue of Hi is zero. Similarly, 
we define H 2 = H 2 — <^i > 0 by offsetting the system Hamiltonian by a constant 
operator Wi. Let pt denote the system state at time t and let pg be the initial state 
of the system at t = 0. We always assume that pg is the ground state of Hi, and so 
we have (Hi)pg = 0. 

The measure of adiabaticity is proposed as follows 
Definition 1 The distance between the final state and the ground state of H 2 is 
measured by 

e={H2)p^. (2) 

Obviously, if the evolution is adiabatic, i.e., pr is the ground state of H 2 , then we 
have e = 0. In particular, e is closely related to the fidelity of the final state and 
ground state in the Schrodinger picture (See Appendix C). A small error e implies 
a large fidelity. 

In this paper we also call e the adiabatic approximation error, as e reflects how 
well we can approximate the evolution as a perfect adiabatic process. 

In this paper we only consider A such that pt, t € [0,T] can be expanded using 
Magnus series in the interaction picture. For more details about the expansion in 
the interaction picture, please refer to Appendix A. If the series expansion is valid 
in the interaction picture, we can transform back to the Schrodinger picture and 
write the evolution of the state as (see Appendix A) 

Pt = e~'^^^\pg-\-R{t)-\-i[pg,X [ dt f{t (3) 

Jo 

where we have ||i?(t)|| = O(A^). A sufficient condition for the Magnus series to 
converge is given by (see Appendix A) 


||AH||/Q^/(t)dt' 


( 4 ) 
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Our aim is to estimate an asymptotic behaviour of e provided A ^ 0. Furthermore, 
we will use the obtained estimate to analyze several cases of the adiabatic theorem 
including those where some difficulties with adiabatic approximation have been 
encountered. 

3 Adiabatic approximation under linear interpolation of the 
Hamiltonian 

The Heisenberg evolution of the expectation of an observable is written as 

(5) 

where H is the system Hamiltonian. Recall that pg = |1)(1|. Since Hi\l) = 

{X{t))p^ is a constant of motion under the action of Hi: 

Hi])p^ = 0 = {-i[X{t),Hi]}p^ (6) 

for any Hermitian operator X(t). 

We will need to study the dynamics of {H 2 )= (i?2(0)ps order to solve 
for e. The time evolution of (H 2 )pt is determined by its generator ^(i? 2 )pt = 

(—i[i?2,^^^(0])pf For linear interpolating function f{t) = y, integration of ^(i? 2 )pt 
over [0,T] results in the following expression (See details in Appendix B): 

m 

( 7 ) 


( 8 ) 

As we noted before, {H2)pt is exactly zero if pr is the ground state of H2- If Pt is 
not the ground state of i? 2 , we can determine the bound on e = {H2)pt from the 
following equality 

{H2)p^-{H2)p, = r{-i[H2,H{t)])p,dt 

Jo 

= / {-i[H2,H{t)])p,dt={H2)p^-{H2)p^. (9) 

Since 


{H2)p^ - {H2)p^ = / {{-i[H2.H{t)])p,)dt 


= f dt[-2{l - - uJi){l\H 2 \i){i\H 2 \l) [ dt f{t ) cos{{uJi - uJi){t - 

Jo Jo 

+ [ dtTr{-iei^i‘[iJ2,(l-/(t))i?i]e-i^^‘i?(t)}. 

1 

A 




2/1 

pT 


t 


dtTr{-iei"i‘[il2, (1 - -)Hi]e-^^-^R{t)} 


H2 — H2 — UJi, 


( 10 ) 
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The error e can be expressed as 


( 11 ) 


With the aid of (8), we can investigate the rate of convergence of e to zero as A 
tends to zero in the case where f{t) defines a linear interpolation, as summarized 
in the following proposition: 

Proposition 1 Assume A 2 > 0 (the ground state of Hi is non-degenerate) and 
suppose f{t) = t/T, which corresponds to the linear interpolation of the Hamilto¬ 
nian. The estimation of e is given by ^ 0{X^), which 

is of the order 0(^4]^) as A —>■ 0. 

Proof 1 Referring to (11) and (9), we need to compute the difference between (8) 
and uji — (l|_ff 2 |l)- First we write (8) as 


(7J2)p. - = E(- 


1 4sin2(A,T/2) 


rT 


A. 




){i\H^\i)m2\i) 


( 12 ) 


dtTr{-ie*^^‘[ff2, (1 - /(t))Hi]e-*^i‘i?(t)} 


-E 

pT 


X^\{l\AH\z)\^ 


A. T^X) 

dtTr{-ze*^^*[ff2, (1 - /(t))Hi]e-*^i‘i?(t)}, (13) 


4A2 sin^(A*r/2)|(l|Ai/|t)p 


by noting that 


E' 

i^i 


X, 


_ l(l|-f^i + XAH — A^|(l|Ai/|i)P 

Xi 


1^1 


i^l 


x^ 


(14) 


Moreover, by the definition of the notation O(-) in Section 2 we can write 

^ 4A" sin=^(A.T/2)|(l|Agp)|" - Ol 

fnf 

Denote H = maxj(t)g(op) ||iJ(t)||. Since 

pT . . rp\ 

II dtTr{-*e*^^‘[i?2,(l-/(t))i^i]e-*^^‘i?(t)}||<—ij2||i?(t)|| (15) 

is 0(A3), we can further write (13) as 


{F[2)pt - {H2)pg — - T 


A^I(1|AHK)|^ 

A, 


-E 

i/1 


4A2 sin2(A,r/2)|(l|Ai/|i)|2 
T2A3 


+ 0(A3 


(16) 


Next we will calculate ui 


1 


— (l|i? 2 |l)- yVe have 


(l|i?2|l) = (l|i?i + XAH\l) = + X{1\AH\1). 


(17) 
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The smallest eigenvalue of H 2 can be calculated using the first-order time- 
independent perturbation theory for non-degenerate system. Assume Hi is the un¬ 
perturbed Hamiltonian and the perturbation is XAH, then the lowest eigenvalue of 
the perturbed Hamiltonian Hi + XAH can be written as series in terms of X and wi 
[28]: 


= a;i + A(l|Ai7|l) - ^ + 0(A3). 


(18) 


Thus we conclude 






(19) 


Comparing (16) and (19), the terms |(i|A^g|a)| gg ^^g g^j^gr e 

is estimated by 


e 


^4A2sin2(A.r/2)|(l|AiJ|*)p , ^,^3, 
^ T^Af 


i#i 


o( 


T^Xl 


), A —>■ 0. 


( 20 ) 


4 Error Estimation for Arbitrary Interpolations 

The approach derived in the previous section can be easily generalized for arbitrary 
given continuous interpolating functions. The generalization can simply be done by 
replacing the linear interpolation function with the given continuous function f{t) 
and then recalculating the double integration 

Ai{T) = -2f dt f dt {1 - ^)Xtf{t) cos{Xfit -t)) (21) 

Jo Jo ^ 

in (7). The error estimation can easily be obtained from the proof of Proposition 1: 


Proposition 2 For an arbitrarily given fit), the error estimation is given by 

e = A^ ^ Al,(r)|(l|AiJ|i)|2 + A^ ^ *^^*^^*'^*" + 0(A3) (22) 

iiii 1^1 * 

as X —y 0. 

Proof 2 e is still calculated by (11), using (i72)pT ~ {H 2 )pg and uji — {H 2 )pg. We 
have 

{H2)p^ - {H2)pg =Y.AfiT)Xfi{l\AH\i)\^ + 0(A3) (23) 


and 


Wi - {H2)pg 


-vv!<1^ + o(a3). 


(24) 
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It must be pointed out that A{T) is very easy to calculate with the aid of any 
softwares that can perform symbolic integration, and therefore it is straightforward 
to apply Proposition 2 to find the error estimation for a given interpolating function, 
as we are going to do in the next section. 


5 Examples 

5.1 Linear Interpolation: f{t) = t/T 

By Proposition 1, the error estimation is e = ^ |(11A|z)pA^+C>(A^) 

as A —^ 0. Since sin^(AiT/2) and AH are bounded, this error term is primarily deter¬ 
mined by y which is the average speed of the variation of the system Hamiltonian, 
and ^ which is the inverse of the energy gap between the ground and i-th eigen¬ 
states of Hi, as A —>• 0. In particular, we have 


lim —r— 

A-o ay 


E 

J/1 


4sin2(A,T/2) 

A? 


|(l|AH|z)|^ 


(25) 


Therefore, when the inverse of the energy gaps are fixed values, the approxima¬ 
tion error e is estimated to be proportional to the square of the average speed of 
the variation of the Hamiltonian, which is (^)^, as A —?> 0. 


5.2 Quadratic Interpolation: f{t) = jT"^ 

Replace /(t) with a nonlinear function /(f) = ^ in (7) and we recalculate the 
integral to be 




1 16sin2(^)-pdT^Af - 8rAiSin(TA0 


+ 


i/i 


i#i 




)|(l|AH|*)|^ 


(26) 


By Proposition 2, for sufficiently small A, the error is estimated to be of order of 
A^: 


^quad 




16sin^(^) + 4r2A2 




(K 2 ^ R6sm^(^) ^ 

iyl y2A5 + xf 


EEEE!l|(i|AH|*)|2 + 0(A3) 

^^|^]|(1|AH|z)P + 0(A3). 

(27) 


That is, in contrast to the linear interpolation case, we have 


lim 




A^O (A)2 


i/1 


A? 


TXj 


(28) 


This calculation shows that if the evolution speed is infinitely slow, then the system 
dynamics is adiabatic during t € [0,T]. However, the scaling of Cguad with respect 
of the square of the average evolution speed y is not as simple as in the linear 
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case, where the scaling of e with respect of is primarily determined by the 
inverse of the energy gaps as A —>■ 0. In the quadratic case, this scaling is primarily 
determined by a complex factor ^ ^ ~ which depends mainly 

on the inverse of the energy gaps {A^} and the inverse of the evolution time T. 

5.3 Interpolation with Decaying Resonant Terms 

Here we assume a linear interpolating function with an additional oscillating term 
that gradually decays to zero. That is, 

fit) = ^ + 5(1 - ^)sin(Act), 

where Ac is the oscillating frequency of the perturbation. Ortigoso observed in [25] 
the inconsistency in the applicability of the adiabatic theorem when the Hamiltonian 
contains resonant terms whose amplitudes go asymptotically to zero. 

Replace /(<) with f{t) = + + g{l — y)sin(Act) in (7) and we recalculate the 
integral to be 


i#l ^ ^ * 


Qi(g,r, Aq Ac) _ 

8A9A2 + 12AjA4-8A5A6 + 2AfA8)- 


(29) 


Qi is a function of four parameters. In particular, we note that each term in (29) 
is well defined for all Ac, including Ac = Aq since as Ac —>■ \i, the i-th term in (29) 
approaches 


[—128sin^(—^) — I6gsin(TAi) + 8gsin(2TAi) + 

+ g‘^{2sm^{TXi) - I) - 2,2T'^\^ - l&gTX^ + 2g‘^T^Xt 

TA 

+ IhgT^A^ sin(TAi) - 16grA*(2 sin^)^) - I) 

- 2g‘^T^XU2sm‘^{TX,) - I) - 2g‘^TXiSm{2TXi)]/32T^Xl 

= - V + ^ sin(rA0 - ^(2 sm\TX,) - I) + Q(T), (30) 

Xi 15 2Ai 15Ai 

where Q(T) is a complicated fraction with T being in its denominator. The error 
resulting from the i-th term is given by 


= l(l|Ai7|*)| 

+ 0(A3) 


2,g^T^X, , ffr2sin(TA,) g^T^{2sm\TX,)-l) , ^2 


16 


2A. 


16A* 


T^Qinij;)" 

(31) 


as A —>■ 0. We have 


A^O (A)2 16 2Xi 16Xi 

( 32 ) 
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The scaling of with respect of is additionally determined by and T^, as 
compared to the quadratic case. This is where adiabatic approximation error may 
not be small if the average evolution speed is slow. In particular by (32), if one 
chooses a comparably large value for T in an adiabatic evolution experiment, the 
adiabatic approximation error may not decrease as expected when one applies a 
slow evolution speed y. 

In order to further illustrate this point, we can heuristically compare the speed of 
convergence of e to zero observed in this case and the quadratic case, as the speed 
of the adiabatic process (A/T) reduces and the evolution horizon T increases. The 
difference in the speed of convergence can be clearly seen using the ratio 

lim I lim — ) = oo. (33) 

T—^oo y(A/T)^0 Cquad J 

Therefore, the rate of convergence considered in this subsection is slower than that 
in the quadratic or linear case, i.e., e goes to zero as A 0 at a much slower rate 
than in the linear interpolation case or the quadratic interpolation case if T is large. 
Furthermore, the larger T is, the slower the convergence. 

6 Conclusion 

In this paper we provide a rigorous analysis of the time-dependent evolution of 
Hamiltonian-varying quantum systems. As we calculated, the adiabatic approxima¬ 
tion error is not proportional to the average speed of the variation of the system 
Hamiltonian and the inverse of the energy gaps in many cases. The results in this 
paper may provide guidelines when applying complicated interpolation for adiabatic 
evolution. 

Appendix A 

The Magnus expansion is proposed to solve the following time-dependent equation 

[29| 


^ = Amity 

The solution of the above equation can be written as 


(34) 


Y{t) = exp(^Ofc(t))y(0), 


(35) 


where the first three terms in the Magnus series {12^, fc = 1, 2, • • •, oo} are calculated 
by 

Hi(t) = f A{ti)dti, 

Jo 

^2{t) = — f dti f dt2[A{ti), A{t2)], 

^ Jo Jo 

^3(i) = ^ / *1 / ^^2 f dt3{[A{ti),[A{t2),A{t3)]] + [A(t3),[A{t2),A{ti)]]). 

0 Jq Jo Jo 

(36) 
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The rest terms in the Magnus series can also be written as the integrals of nested 
commutators. 

The dynamical equation of the quantum state in interaction picture is given by 

(37) 

where |^/’(t)) = Applying the Magnus expansion to (37) yields 

|V’/(t)) = (1 - iA / /(t')Ai7e“^^^* + i?o(t))|^/’(0)), (38) 

Jo 

where Ro(t) includes all the higher-order terms as determined by {flfc}. Obviously, 
||i?o(t)|| is of the order O(A^). Transforming back to Schrddinger picture we obtain 
the expression for the density operator as 

Jo 

+ Ro{t))pg{l + iX [ dt f{t + Rl{t))e^^^* 

Jo 

= e-i^^‘(pg + i7(t) + i[p3,A / dt'ei"i‘'/(t')A77e-i"i‘'])ei"i*. (39) 

Jo 

Obviously, ||i?(t)|| is also of the order O(X^). 

An explicit condition for the Magnus series to converge is given by [29] 

||Ai7|| rA/(t)dt <^. (40) 

Jo 


Appendix B 

The derivative of (772 )pt is calculated as 


-b 


-b 


Hi + /(t)(ff2 - Hi)])„ = {-ilHi.Hi + m{H2 - Hi)]),, 

-([ei^^‘[i72,(l-/(t))i7i]e-i^^‘,A [ dt f{t )AHe-'^^^^']), 

Jo 

Tr{-iei^i‘[i72, (1 - R{t)} 

-([ei^^‘[i72,(l -/(t))i7i]e-i^^‘, f dt f{t ]) 

Jo 

Tr{-iei^i‘[i72, (1 - R{t)}, 


/Pg 


Pg 


(41) 



where we made use of the relation (6) and XAH = H 2 — Hi -\- uii — uii. Calculating 
(41) further leads to 


— lUJlt 


= -(1 - dt f{t)H2e-'^ 

Jo 

+ Tr{-iei^i‘[iJ2, (1 - 

= -(l-/(t))(ei‘"i‘iJ2-ffi V|z)(z|e-i^^* f dt f{t)H2e-'^^^* 

Jo 

Jo 

+ Tr{-iei^i‘[i?2, (1 - f{t))Hi]e-'^^^^R{t)} 

= -(l-/(t))^(w.-ccl)(i^2|^)(^|i^2)p, f dt 

z/l -^0 

+ Tr{-iei^i‘[iJ2, (1 - fit))Hi]e-'^^^^R{t)} 

= -2(1 - /(t)) V(Wi - Wi)(l|i? 2 N)(*|-ff 2 |l) [ dt f{t)cos{{uJi-uJi){t -t)) 

.-^1 Jo 


i^l 


+ Tr{-iei^i‘[iJ2, (1 - /(t))i?i]e-i^^*i?(<)}. 

With the linear interpolating function f{t) = the direct integration of (42) over 
[0,T] gives (8). 


Appendix C 

The state of the system will remain a pure state during the evolution. Therefore, 
we can express the final state as px = |V’)(^I with W = )• Using this 

expression, the error measure e defined by (2) can be written as 

N N 

e = {H2)pr = {m2\i’) = E ^ |c.Pa;. (43) 

i—2 i—2 

The fidelity of the final state and the ground state |1 ) is calculated by 


nw),ii')) = Vkv'|i')f = %/rf= 


\ 


1 - 


N 

El 

2 = 2 




(44) 
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